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Abstract 

Arrays of identical limit-cycle oscillators have been used to model a wide variety of pattern- 
forming systems, such as neural networks, convecting fluids, laser arrays, and coupled biochemical 
oscillators. These systems are known to exhibit rich collective behavior, from synchrony and 
traveling waves to spatiotemporal chaos and incoherence. Recently, Kuramoto and his colleagues 
reported a strange new mode of organization — here called the chimera state — in which coherence 
and incoherence exist side by side in the same system of oscillators. Such states have never 
been seen in systems with either local or global coupling; they are apparently peculiar to the 
intermediate case of nonlocal coupling. Here we give an exact solution for the chimera state, for 
a one-dimensional ring of phase oscillators coupled nonlocally by a cosine kernel. The analysis 
reveals that the chimera is born in a continuous bifurcation from a spatially modulated drift state, 
and dies in a saddle-node collision with an unstable version of itself. 
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I. INTRODUCTION 



The chimera state 



A fascinating; spatiotemporal pattern was reported recently by Kuramoto, Battogtokh and 

Ho, 



Shima 



11 



18| . While studying arrays of identical limit-cycle oscillators that are coupled 
nonlocally, they found that for certain choices of parameters and initial conditions, the array 
would split into two domains: one composed of coherent, phase-locked oscillators, coexisting 
with another composed of incoherent, drifting oscillators. The coexistence of locking and 
drift was robust. It occurred in both one and two spatial dimensions, and for various kinds 
of oscillators, including the Fitzhugh-Nagumo model, complex Ginzburg-Landau equations, 
phase oscillators, and an idealized model of biochemical oscillators. 

It's important to appreciate how unexpected this coexistence state was. Nothing like it 
had ever been seen before, at least not in an array of identical oscillators. Normally, identical 
oscillators settle into one ofafew basic p*™ flfl0. The s impl esUs s^c h „ h 
all oscillators moving in unison, executing identical motions at all times. Another common 
pattern is wave propagation, typically in the form of solitary waves in one dimension, spiral 
waves in two dimensions, and scroll waves in three dimensions. The common feature in 
these cases is that all the oscillators are locked in frequency, with a fixed phase difference 
between them. At the opposite end of the spectrum is incoherence, where the phases of 
all the oscillators drift quasiperiodically with respect to each other, and the system shows 
no spatial structure whatsoever. And finally, one sometimes sees more complex patterns, 
including modulated structures, spatiotemporal chaos, intermittency and so on. 

What was so odd about the coexistence state is that two of these patterns (locking and 
incoherence) were present in the same system, simultaneously. This combination of states 
couldn't be ascribed to the simplest mechanism of pattern formation — a supercritical insta- 
bility of the spatially uniform oscillation — because it can occur even if the uniform state 
is linearly stable, as indeed it was for the parameter values used by Kuramoto and his 
colleagues. Furthermore, it has nothing to do with the classic partially locked/partially 
incoherent states that occur in populations of non-identical oscillators with distributed nat- 
ural frequencies |3| . There, the splitting of the population stems from the inhomogeneity 
of the oscillators themselves; the desynchronized oscillators are the intrinsically fastest and 
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slowest ones in the tails of the distribution. In contrast, for the system studied by Kuramoto 
et al., there is no distribution of frequencies. All the oscillators are the same, and yet they 
still break apart into two groups of utterly different character. 

Because the coexistence state involves two seemingly incompatible forms of dynamical 
behavior, we will henceforth refer to it as "the chimera state," inspired by the mythological 
creature composed of a lion's head, a goat's body, and a serpent's tail. Today the word 
chimera is used more generally to indicate something made up of incongruous parts, or 
something that seems wildly improbable or fantastical. 

Figure^shows a realization of the chimera state in the simplest setting, a one-dimensional 
ring of phase oscillators U|- The governing equation is 



dt 



u- I G(x — x') sin [<p(x, t) - <f)(x', t) + a] dx' (1) 
Jo 



where <p(x, t) is the phase of the oscillator at position x at time t. The space variable x runs 
from to 1 with periodic boundary conditions, and should be regarded as an angle on a circle 
(mod 1). The frequency uj plays no role in the dynamics, in the sense that one can set uj = 
without loss of generality by redefining 0^0 + ut, without otherwise changing the form of 
equation ((TJ). The kernel G(x — x') provides nonlocal coupling between the oscillators. It is 
assumed to be even, non-negative, decreasing with the separation \x — x'\ along the ring, and 
normalized to have unit integral. Specifically, Kuramoto and Battogtokh 1 1 11 ] used an 
exponential kernel G(x — x') oc exp (— k\x — x'\). Then, for parameter values a = 1.457 and 
k = 4 and suitable initial conditions (to be discussed in detail in Section ITTT]). the system 
evolves to the chimera state shown in Fig. ^ 

In this snapshot of the instantaneous phases, two distinct regions are conspicuous. The 
oscillators near x = (and equivalently, x — 1) are phase-locked. All of them move with 
the same constant frequency; in a frame rotating at this frequency, they would all look 
frozen. The smoothness and flatness of the graph of <f>(x) in this region indicates that these 
oscillators are coherent as well, i.e., they are nearly in phase. 

Meanwhile, the scattered oscillators in the middle of Fig. ^are drifting, both with respect 
to each other and with respect to the locked oscillators. Their motion is strongly nonuniform. 
They slow down when they pass near the locked pack — which is why the dots appear more 
densely clumped at this phase — and then speed up as they lap it. 
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FIG. 1: Phase pattern for a typical chimera state. Here k = 4.0, a = 1.45, N = 256 oscillators. 
Equation Q was integrated with fixed time step dt = 0.025 for 8,000 iterations, starting from 
4>(x) = 6exp — 30(x — i) 2 r(x), where r is a uniform random variable on 
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FIG. 2: Phase pattern for a typical chimera state shown on the torus. Azimuthal angle indicates 
spatial position x. Phase <j) is constant along lines of latitude; the outer equator of the torus 
corresponds to (ft = 0. 



B. Puzzles 

When we first learned about the chimera state by reading Ref. Q|, we were amazed by 
it. How could such a thing even be possible? 

In fact, a little thought showed that it was provably impossible in two special cases that 
had been studied previously: 

• Global coupling: Chimera states can't occur for Eq. (0) with G(x) = 1 and any choice 
of < a < 7r/2, because a Lyapunov function exists for this case, demonstrating that 



almost all solutions are attracted to the in-phase oscillation 23|, |24 ] . 



• Sine coupling: If a = 0, corresponding to a pure sine coupling in Eq. (JTJ), chimera 
states are impossible for any even kernel G of any range. This follows because Eq. (JIJ 
becomes a gradient system in the frame rotating at frequency us. Hence all attractors 
must be fixed points, corresponding to phase-locked solutions in the original frame, 
thus ruling out the possibility of coexisting drift. 

So the coexistence phenomenon must somehow rely on a conspiracy between a^O and the 
non-global nature of the coupling. But how, exactly? 

And for that matter, is the chimera state born as soon as a 7^ 0, or at some value of a 
bounded away from zero? In dynamical simulations like that shown in Fig. ^ stable chimera 
states are observed only when a is close to, but slightly less than, tt/2. Does that mean that 
these states don't exist for smaller a, or is just that their basins of attraction shrink as a 
decreases? 

Furthermore, what is the genealogy of the chimera state, in the sense of bifurcation 
theory? Is it born out of the vacuum, as a pair of stable and unstable versions of itself? 
Or does it emerge when some more familiar attractor loses stability? For instance, does 
it bifurcate off the fully incoherent state, in which oscillators are uniformly scattered and 
drifting around the circle at every x? That seems unlikely, since the phase pattern shown 
in Fig. El looks pretty far from total incoherence; even its drifting oscillators show some 
clumping in phase. So maybe the chimera state branches off the uniform in-phase state? 
But how can it, given that the in-phase state is linearly stable for all \a\ < |? 

Motivated by these puzzles, we have tried to understand where the chimera state comes 
from and to pinpoint the conditions that allow it to exist. A brief report of our findings 
appeared in 



C. Broader significance 

Aside from the questions it raises, we believe the chimera state is also more broadly 
significant for nonlinear science, for two reasons. 

First, it exemplifies the surprises that lurk in nonlocally coupled systems. As Kuramoto 



and his colleagues have pointed out 



tnat lurk m nonlocally 



, nonlocal coupling is a relatively 
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dark corner of nonlinear science in general, and nonlinear oscillator theory in particular. 
Most previous work on coupled oscillators has focused on local coupling, where the inter- 
actions are assumed to be solely between nearest neighbors, or global coupling, where each 
oscillator interacts equally strongly with all the others. The intermediate case of nonlocal 
coupling is natural to explore next, and has already revealed some interesting new forms of 
dynamical behavior P. 

From a more applied perspective, nonlocal coupling is important to investigate because 
it arises in diverse systems throughout physics, chemistry, and biology. Examples include 
Josephson junction arrays [l^, chemical oscillators epidemiological models of disease 
spread |3| , and the neural networks underlying the patterns on mollusc shells ((J Q , lo- 
calized neural "bump" states 
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Second, the chimera state is by no means an oddity restricted to Eq. (JI}. On the contrary, 
it was first seen in simulations of the complex Ginzburg-Landau equation with nonlocal 
coupling , a fundamental model in the study of pattern formation. That equation in 

turn can be systematically derived from a wide class of reaction-diffusion equations, under 
particular assumptions on the local kinetics and diffusion strength that render the effective 
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coupling nonlocal 



181 . 1221 ] . Under an additional assumption that the coupling is 
also sufficiently weak (in a precise sense), Shima and Kuramoto [l8[ show that the original 
react ion- diffusion system can be further reduced to a phase equation of the universal form 



Of 



J dr'G (r 



u - / dr'G (r - r') sin \6(r) - 6(r') + a] 



where r labels the position of the oscillators and the kernel G decays exponentially with 
distance: G{r — r') oc exp (— n\r — r'\) . But this is just Eq. ((TJ), if the space is one- 
dimensional. So there is good reason to expect that the coexistence phenomenon will have 
some generality. 

For example, in two dimensions, the coexistence of locked and drifting oscillators mani- 
fests itself as an unprecedented kind of spiral wave: one without a phase singularity at its 
center [l^, Q . Instead, the oscillators in the core are found to be completely desynchro- 
nized from each other and from the uniform rotation of the spiral arms. In effect, the core 
oscillators mimic a phase singularity by being incoherent. A better understanding of the 
one-dimensional case might shed light on this remarkable new form of pattern formation. 
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II. SUMMARY OF PRIOR RESULTS 



We begin by reviewing the results of Kuramoto and Battogtokh |l(J,|ll||. After uncovering 
the chimera state in their simulations of Eq. (jl)). they were able to explain much of its struc- 
ture analytically Their elegant approach is a generalization of Kuramoto's self-consistency 
argument for globally coupled oscillators |i| ^| . 

In this approach, one first transforms by seeking a rotating reference frame in which 
the dynamics become as simple as possible. Let Q denote the angular frequency of this 
rotating frame (to be determined later, in the course of solving the problem), and let 

9 = - ttt 

denote the phase of an oscillator relative to this frame. Next, introduce a complex order 
parameter Re te that depends on space and time: 

R(x, t)e i6 ^ = f 1 G(x - x') e^'^dx' . (2) 

J 

To see what this order parameter means intuitively, note that the integral on the right 
hand side of (j2J) performs a running average of e td over a window centered at x, with a width 
determined by the width of the kernel G. Thus < R(x, r) < 1 can be viewed as a measure 
of the local phase coherence at x, and 9(x, t) represents the local average phase. These 
two average quantities provide macroscopic proxies for the overall state of the continuum of 
oscillators. 

The real virtue of introducing the order parameter, however, is that we can now rewrite 
the governing equation (JTJ) as 

39 

— = LU-n-Rsm[9-@ + a] , (3) 

which makes it look as if the oscillators have decoupled, though of course they are still 
interacting through R and G, to which they each contribute through (J2J). This observation 
suggests that the problem can be attacked by the self-consistency arguments of mean-field 
theory, even though it is not globally coupled. 

Now comes the key step. Suppose we restrict attention to stationary solutions, in which 
R and depend on space but not on time. Now the equations truly do decouple, in the 
following sense. One can easily solve for the motion of the oscillator at each x, subject to 
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the assumed time-independent values of R(x) and Q(x). The oscillators with R(x) > — 
asymptotically approach a stable fixed point 9*, defined implicitly by 

io -Q = R(x) sin [6* - &(x) + a] (4) 

The fact that they approach a fixed point in the rotating frame implies that they are 
phase-locked at frequency Q in the original frame. On the other hand, the oscillators with 
R(x) < \uj — fi| drift around the phase circle monotonically. To be consistent with the 
assumed stationarity of the solution, these oscillators must distribute themselves according 
to an invariant probability density p(0). (To ease the notation here and elsewhere, we have 
suppressed the dependence on x whenever it's clear from context.) And for the density to 
be invariant, the probability of finding an oscillator near a given value of 9 must be inversely 
proportional to the velocity there. From (J3J), this condition becomes 

^J(u - Q) 2 - R 2 

P ^ = 2n\cu-tt- Rsm(9-e + a)\ ^ 

where the normalization constant has been chosen such that J^ n p(9)d9 = 1. 

The resulting motions of both the locked and drifting oscillators must be consistent with 
the assumed time-independent values for R(x) and Q(x). To calculate the contribution that 
the locked oscillators make to the order parameter (J2J), observe that 

u — Q 



R 



sin(fl* -Q + a) 

JR 2 -(cu- n) 2 

cos(#* - 6 + a) = ±-¥ (6) 

R 

for any fixed point of (j3J). One can check that the stable fixed point of corresponds to 
the plus sign in (0). Hence 



IR 2 - (Lu-n) 2 + i(uj-n) 
exp [i(9* -© + «)] = ^ i ^ i (7) 

which implies that the locked oscillators contribute 

^R 2 -(u- Q) 2 + i(u - Q) 

R 

(8) 

to the order parameter (J2J. Here the integral is taken over the portion of the domain where 
R(x') > \u-Q\. 
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J dx'G{x - x') exp [i9*(x')} = e~ ia J dx'G{x - x') exp \iQ{x') 



Next, to calculate the contribution from the drifting oscillators, Kuramoto and Battog- 
tokh jlOj, replace exp [i6(x')] in (J2J with its statistical average J^ 7T exp(i6)p(8)d8. Using 
(J5J) and contour integration, they obtain 

J* exp{i6)p{6)d6 = -i (u - Q - J(u - Q) 2 - R 2 ^ . 

Therefore the contribution of the drifting oscillators to the order parameter is 

r r ■ r u — Q — J(u - n) 2 - R 2 (x') 
/ dx'G{x-x') / exp(i9)p(6)d9 = ie~ ia / dx'G(x-x') exp [i®(x')} v — 

where now the integral is over the complementary portion of the domain where R(x') < 
\uj — 

Notice something curious: the integrand on the right hand side of the drifting contribution 
is exactly the same as that found earlier in dHl for the locked contribution: only their domains 

□n 

differ. (This coincidence is not mentioned in [10,lU(.) To see that the two expressions agree, 
note that 



r 2 - (w - n) 2 + i(u-n) = i\ u-n-^/(u- n) 2 - r 2 

as long as we choose the branch corresponding to the square root of a negative number. 
Hence the two contributions can be combined into a single integral, yielding a slightly 



more compact version of the self-consistency equation derived in 
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. r i iv-Q- J(LU-tt) 2 -R 2 (x') 

R(x) exp \iG(x)] = ie~ ia / G(x - x') exp \iQ(x')\ V — ' — dx' . 

Jo R(x') 

To ease the notation a bit more, let 

3 = a 

H 2 

A = a; - Q . (9) 
Then the self-consistency equation becomes 



■ a r 1 - a , „A-JA 2 -R 2 (x') 

R ( x )e lQ W = e lf3 / G(x - x')e l ^ x ) ^— —dx' . (10) 

Jo R(x') 

Equation ()10|) is to be solved for three unknowns — the real-valued functions R(x) and 
Q(x) and the real number A — in terms of the assumed choices of f3 and the kernel G(x). 
Notice that although u itself is arbitrary up to a constant, and hence so is Q, their difference 
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to — is physically meaningful; it is determined by the condition that the long-term dynamics 
become stationary in the frame rotating at frequency Q. 

Initially, we couldn't see how to solve the self-consistency equation (|l(Jj) numerically. 
We wrote to Kuramoto for advice, and he described an iterative scheme to determine the 
functions R(x) and G(x), based on initial guesses obtained from the dynamical simulations. 
The idea behind the scheme is that the current estimates of R(x) and G(x) can be entered 
into the right hand side of (jlOJ) . and used to generate the new estimates appearing on the 
left hand side. 

Still, that leaves open the question of how to determine A. We seem to have only 
two equations (given by the real and imaginary parts of Eq. (jl(Jj) ) for three unknowns. 
Fortunately, a third equation can be imposed to close the system. Because (ITU]) is left 
unchanged by any rigid rotation Q(x) —>■ Q(x) + 0o, we can specify the value of O(x) at any 
point x we like; this freedom is tantamount to choosing an origin in the rotating frame. A 
natural choice would be to demand 0(0) = 0, but as we'll see in Section \W\ another choice 
turns out to be more convenient. 

Kuramoto and Battogtokh ^j, 11] confirm that the self-consistency approach works: their 



results from numerical integration of the dynamical equations (0) match those obtained by 
solving the self-consistency equation (fTUj) iteratively. 

Figure El shows the chimera state along with the graphs of R(x) and Q(x) for the param- 
eters used in Fig. 21 The curves in Fig. E^b) and|Hfc) are periodic and reflection-symmetric. 
In fact, they resemble cosine waves, which made us wonder whether (jlO)) might have a simple 
closed-form solution, perhaps in some perturbative limit as a parameter tends to zero. To 
see where such a limit might come into play, we hoped to first replicate the simulations of 



Kuramoto and Battogtokh 
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111 ] and then to explore parameter space more widely. 



III. A FIRST ROUND OF SIMULATIONS 

Unfortunately, we couldn't find the chimera state in our early simulations of Eq. (0). No 

matter how we started the system, it always converged to the in-phase state. In the report 

I I 

that announced the chimera state, Kuramoto |H| does not give precise details of the initial 



condition he used. He describes it as a "suitable single-humped initial phase pattern" 
p. 219] which we incorrectly took to mean something like <p(x, 0) = a + b cos x or e~ acoE 
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(a) 




FIG. 3: Chimera state and order parameter curves for the exponential kernel G(x) oc exp (— k|x|), 
as used by Kuramoto and Battogtokh 2 5 11 1- Parameters are the same as those in Fig. ^ (a) 
Phase pattern for chimera state, (b) Local phase coherence R(x), computed from (J2J). Locked 
oscillators satisfy R(x) > A. (c) Local average phase Q(x), computed from (J2J). 

Eventually, we asked Kuramoto for help (again!), and he kindly explained what he meant. 
(He also sent us his paper with Battogtokh, where the description of the initial condition 
is more explicit.) At each x, a uniform random number <p(x, 0) is chosen within some 
interval whose width varies with x in a single-humped fashion. Specifically, the width is 
narrowest near x = (mod 1), meaning that the oscillators are most nearly in phase there, 
initially. As x increases toward the diametrically opposite point of the domain at x = |, the 
phases are scattered progressively over larger and larger regions on the phase circle (meaning 
the oscillators are placed more and more incoherently there, initially). The effect of this 
procedure is to give the system a jump-start, by placing it in a partially coherent/partially 
incoherent state to begin with. 

To be more precise, Kuramoto used a random distribution with a Gaussian envelope: 

r / , \ 2" 

(p(x, 0) = 6 exp 



r(x), where r(x) is a uniform random number on the interval 
— \ < r < \ ■ For the parameters used in Fig. this initial condition indeed evolves to the 
chimera state reported in 

Then we ran simulations to see how far this state could be continued by decreasing cv, 
knowing that it would have to disappear or lose stability somewhere before a = 0. To track 
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its fate along the way, we also computed several statistics: 

1. the spatial average of R(x), given by (R) = Jo 1 R{x)dx- 

2. the amplitude of R(x), defined as R amp = R ma x — Rmin\ 

3. /drift, the fraction of the spatial domain occupied by drifting oscillators; 

4. the difference A = uo — Q between the nominal frequency of individual oscillators and 
their collective frequency when locked; and 

5- A max , the largest value of the time-averaged drift velocity relative to the rotating 
frame. This quantity measures the average speed of the fastest drifting oscillator. From 
Q, it can be calculated as maxj. |A(a;)| = max x A 2 — R 2 (x)\, where the maximum 
is taken only over the drifting oscillators. 

Figure |U shows how /,-[ r ift var i es when k is held constant but a is changed smoothly. We 
generated similar graphs for each of the statistics mentioned above, and all showed a jump 
to the uniform synchronized state as a decreased below some critical value a c . From these 
results it appeared that when k = 4.0, the chimera state ceased to exist somewhere around 
a c ~ 1.37. The transition seemed to be discontinuous, which suggested that a c couldn't be 
calculated by a naive perturbation expansion. If it was to be calculable at all, something 
more subtle would be required. 

The next step was to investigate how these results vary with k. Recall that the kernel 
in (JTJ) is G{x) = Cexp(— k\x\), so 1/k sets a characteristic length scale. Roughly speaking, 
it is the distance over which the nonlocal coupling is substantial. So the limit n — > + 
corresponds to global coupling G(x) = 1. This can also be checked directly, noting that the 
normalization constant for the exponential kernel on the circle is given by C = f (1— e _K / 2 ) -1 . 

Figure El shows a rough contour plot of R amp in the (a, k) parameter plane. Crude as 
this plot is, its message is still clear. The stable chimera state evidently lives in a wedge 
in parameter space, bounded on one side by the line a = | and on the other by a curve 
a = a c (n) that is nearly a straight line. By its very shape, the picture directs our attention 
to the corner of the wedge, to the simultaneous limit as a — > | from below and k — > from 
above. Apparently something crucial happens in that corner — the chimera state is born 
there. And so this is where perturbation theory should be conducted. 
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FIG. 4: The fraction of oscillators drifting as the coupling parameter a varies. Here k = 4.0, 
N = 256 oscillators, dt = 0.025 for 100,000 iterations. 

To check that the wedge of Fig. El was not an artifact of the exponential kernel assumed 
above, we also calculated the corresponding contour plots for the cosine kernel 

1 



G(x) = — (1 + Acosx) 
2ir 



(11) 



where < A < 1. Here the spatial domain has been redefined to — 7r < x < tt for 
convenience, and to bring out its ring geometry and the reflection symmetry of the chimera 
state. Figure |H1 confirms that the cosine kernel gives a similar chimera state to that for the 
exponential kernel used above, while Fig. [7| demonstrates that the wedge in parameter space 
is preserved as well. All that is reassuring, because as it happens, the cosine kernel also has 
the pleasant property that it allows the self-consistency equation to be solved analytically. 



IV. AN EXACTLY SOLVABLE CASE 

From now on, let G(x) be given by the cosine kernel (jllj) . and let the spatial domain 
be — 7r < x < 7r with periodic boundary conditions. For this case, we'll show that the 
functional form of the order parameter can be obtained exactly, which in turn yields the 
explicit x-dependence of R(x) and 0(x). All the resulting expressions, however, still contain 
two unknown coefficients, one real and the other complex, which need to be determined self- 
consistently. In this way, the two unknown functions in the self-consistency equation are 
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FIG. 5: Amplitude of the curve R(x), depicted as a contour plot in parameter space, and 
calculated by averaging over the instantaneous R curves during numerical integration. Here 
G{x) oc exp(— N = 80 oscillators, the integration time step is dt = 0.025, and integration 
continued for 20,000 iterations. Lighter colors indicate smaller amplitude; lightest is R am p = 0.0 
and darkest is R am p = 0.13. 

exchanged for two unknown numbers — a drastic reduction in the difficulty of the problem. 
The self-consistency equation ftTHft is 

R(x)e l@{x) = e i/3 r G(x - x')e i6 ^h(x')dx' (12) 

where we've introduced the notation 



A - J A 2 - R 2 (x') 

ft M = V) ' < 13 > 

Let angular brackets denote a spatial average: 

(f) = ^ r f(x f )dx'. 

LlX J — 7T 

Then, substituting the cosine kernel (|11|) into (|12|) and expanding G(x — x') by a trigono- 
metric identity, we find 

Re t0 = — / [l + Acosxcosx' + Asmxsmx']h(x')e ( x >dx' 

2% J-n 

= e if3 (he ie ) + e i/3 A (he ie cosx 1 ) cosx + e^A (he ie sin a/) sinx 

= c + a cosx (14) 

14 
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FIG. 6: Chimera state and corresponding order parameter curves for the cosine kernel, shown in 
the same format as Fig. El and qualitatively similar to it. Parameters are A = 0.995, f3 = 0.18, 
N = 256 oscillators; equation Q was integrated with fixed time step dt = 0.025 for 200,000 
iterations, starting from <j>(x) = 6r exp(— 0.76x 2 ), where r is a uniform random variable on 
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where the coefficients c and a must satisfy their own version of the self-consistency equations, 
now given by 

c = e ip (he 10 ) (15) 

and 



Ae if3 (he ie cos 



x 



(16) 



Note that the coefficient of sin x vanishes in (|14|) . This follows from the assumption that 
R(x') = R(-x') and Q(x') = 0(— x'), as suggested by the simulations; then h(x') in (fTHj) 
is also even, and so the integral (he 10 sin x') in (|T4~jl vanishes by oddness. As we'll show 
next, this assumption of reflection symmetry is self-consistent, in the sense that it implies 
formulas for R(x) and Q(x) that indeed possess this symmetry. 

For example, to calculate R(x) in terms of the unknown coefficients a and c, observe that 



R 2 



(Re ie )(Re- ie ) 

(c + a cos x) (c* + a* cos x) 

|c 2 | + 2Re (ca*) cosx + \a\ 2 cos 2 x 



(17) 



which is an even function, and which also helps to explain why the graph of R(x) in Fig. El 
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FIG. 7: Contour plot of R a mp for the chimera state with cosine kernel. Note the similarity to 
Fig. El for the exponential kernel. Here G{x) = o^(l + Acosx), N = 80 oscillators, the integration 
time step was dt = 0.025, and integration continued for 30,000 iterations. Same color scale as 
Fig. El 

resembled a cosine wave. 

Likewise, 0(x) is an even function reminiscent of a cosine because 

R(x) sin Q(x) 



tan 6 (a;) 

R{x) cos 

Ttyi I n\ -L Tin ( n i r*ric 'r 

(18) 



Re (c) + Re (a) cos x 

Another simplification is that c can be taken to be purely real and non-negative, without 
loss of generality. This follows from the rotational symmetry of the governing equations. 
In particular, the self-consistency equation (J12)) is left unchanged by any rigid rotation 
Q(x) — > Q(x) + 0o- Thus we are free to specify any value of Q(x) at whatever point we like. 
The most convenient choice is to set 

e(!'i-o. 



.2. 

Then at that value of x the equation Re lB = c + a cos x reduces to 



.2. 

Since R is real and non-negative, so is c. Hence, we take 



lm(c)=0 (19) 
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from now on. 

The final step in closing the equations for a and c is to rewrite the averages in (j!5|) and 
(J 16)1 in terms of those variables. To do so, we express he t@ as 

h 



he ie = (Re ie ) 



R 



A- J A 2 -R 2 (x) 
(c + acosx) — 



R 2 (x) 

A - ,/A 2 - R 2 (x) 

= i~i — (20) 

c + a* cos x 

where we have used (j!7)) and the real-valuedness of c to simplify the second line above. 
Inserting (fTTj) and (f2Uj) into (fTB]) and (fTo^) . we obtain the desired self-consistency equations 
for a and c: 

i/9 / A - (A 2 - c 2 - 2Re(a) ccosx - |a| 2 cos 2 a;)2 \ 
\ c + a*cosa: / 

, , fl / A — (A 2 — c 2 — 2Re (a) ccosx — lal 2 cos 2 x) 5 \ . , 

a = Ae lP ( — u — cos x ) . 22 

\ c + a* cos x / 

This pair of complex equations is equivalent to four real equations for the four real unknowns 
c, Re (a), Im (a), and A. The solutions, if they exist, are to be expressed as functions of the 
parameters j3 and A. 



V. CLUES BASED ON NUMERICS 

Before we plunge into the details of solving equations (|2"Tj) and (JUJ) simultaneously, let's 
pause to remember what we're trying to do. 

We want to understand where the chimera state lives in parameter space and how it 
bifurcates. Guided by the simulations of Section UTT] we expect that (|2*Tj) . (|2~2"|) should have 
chimera solutions throughout the wedge-shaped region of parameter space shown in Fig. El 
Assuming that's true, we hope that these solutions will continue all the way down to the 
corner (a, A) = (f,0), corresponding to (/3,A) = (0,0), where might be able to analyze 
them with perturbation theory. 

Our strategy, then, is to start by finding one solution to (J2T|) . (l22|) . by any means possible, 
for parameter values anywhere in the wedge. Having found this solution, we can use it as a 
base point for a numerical continuation method. Then we proceed to dive into the corner, 
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following a straight line through parameter space between the base point and the corner. 
In this way we convert the problem to a one-parameter study of the solutions of (|2ip. (|22j) . 
Sufficiently close to the corner, we expect that the solutions will display some sort of scaling 
behavior with respect to the parameter. That scaling will then suggest clues about the right 
ansatz for a subsequent perturbation calculation. 

So first we have to come up with a chimera solution to (|21jh (|22|). It's not just a matter 
of plugging the equations into a standard root-finding package. The trouble is that these 
equations also have other solutions that we're less interested in, and we don't want the 
numerical root-finding scheme to converge to them instead. 

In particular, the in-phase solution, where all the oscillators are locked at the same phase 
and none of them are drifting, has a large basin of attraction that competes with that of the 
chimera state. To see what values of a, c, and A correspond to the in-phase state, note that 
when 4>(x,t) = (p(x',t) for all x and x 1 , Eq. (JTJ) implies (f>(x,t) = <fi + (u — sina;)£. Hence 
R = 1 and therefore c = 1 and a = 0. And because Q = u — sin a = uj — cos (3, we have 
A = uj — Q = cos (3. Thus 

(a, c, A) in _ phage = (0, 1, cos/3) . (23) 

It's easy to check that this satisfies (|21j). (J22|) for all values of A and (3. 

To reduce the chance that the root-finder will converge onto this in-phase state, we need 
to concoct an initial guess that's very close to a genuine chimera state. To find one, we 
numerically integrated Eq. using the cosine kernel, and fit the resulting graphs of R(x) 
and O(x) to the exact formulas (j!7|) and ()18|) . to estimate the values of a and c. The 
frequency difference A was obtained directly from the simulation, by setting oj = and then 
computing A = uj — Q = — fi, where Q is observable as the collective frequency of the locked 
oscillators. 

In this way we estimated a = 0.156 — 0.072i, c = 0.591, A = 0.720 for the stable chimera 
state at parameter values A = 0.99, j3 = 0.081. We fed this starting guess into the Matlab 
root-finding and numerical continuation program MatCont j^] and found rapid convergence 
to a = 0.162 — 0.051i, c = 0.588, A = 0.723. From there, we could continue the solution in 
either A or (3 or some combination, as we saw fit. 

This approach enabled us to track the chimera state throughout parameter space, until 
it disappeared along a critical curve corresponding to the boundary of the wedge shown 
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FIG. 8: The region of parameter space in which the chimera state exists. Solid line, exact boundary 
determined by numerical solution of 1)21(1 and (|22|): dashed line, leading order approximation to 
this boundary obtained by perturbation theory (see text). 

earlier. The results of this calculation are shown in Fig. |HJ As expected, the boundary of 
the region is nearly a straight line, and it extends down to the origin. 

VI. PERTURBATION THEORY 

The next step is to look for scaling laws to guide our perturbation calculations. Figure 01 
shows the results of numerical continuation starting from (/3, A) = (0.08, 0.99) and moving 
along the line A = 12.375/3 towards the origin, all the while remaining within the wedge 
shown in Fig. |SJ The observed behavior of the variables along that line suggests the following 
ansatz near the origin: 

A ~ 1 + A x e + A 2 e 2 
c ~ 1 + cie + c 2 e 2 
Re (a) ~ ue 2 

Im(a) ~ ve 2 (24) 

where we have introduced e = A as the small parameter. 

Next, we assume that this ansatz continues to hold along other lines through the origin. 
Such lines can be parametrized as 

A = e, 
P = Pie, 
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FIG. 9: Scaling laws near the origin in parameter space, along the line A = 12.375/3. Data were 
collected from numerical continuation of a known chimera state, for an ensemble of parameter 
values. Approximate fits were then determined by least-square regression, (a) Scaling of real and 
imaginary parts of a; (b) Zoom of panel (a) near origin in parameter space. Note that curves are 
quadratic; (c) Linear scaling of real-valued variable c; (d) Scaling of A. Note that A scales linearly 
for small values of A (purple). 

where A and (3 tend to zero simultaneously as e — > 0. Here f3\ is a free parameter inversely 
related to the slope of the lines. Thus the asymptotic shape of the wedge in Fig.|HJ sufficiently 
close to the origin, will be determined from the maximum and minimum values of f3\ for 
which a perturbative solution exists. 

Substituting the ansatz (}24"J) into the self-consistency equation pTf for c, and retaining 
only terms up to O (y/e) gives 

l + 0(e) = (iW ^-P+f*- 1 - 2 ^* ) 

\ 1 + cie / 



= 1- V2\JAi- Cl V^+0{e) , 

implying that 

A x = ci . (25) 

Now we retain terms up to O (e) on both sides, and apply Eq. (|25|) whenever necessary 
to cancel terms. At this order, Eq. ()21j) becomes 



l + cie= 1 + eV2(J(A 2 - c 2 ) -ucosx) . (26) 
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To simplify notation, let 



5 = A 2 - c 2 . 



(27) 



After breaking up the previous expression (J26|) into two equations for the real and imaginary 
parts, and equating terms of O (e) , we get 



Repeating the same expansion to O (e) in the self-consistency equation (|22j) for a yields two 
analogous equations: 



The equations (|2~% )) - (J3*T]) form a closed system for the variables (c±,u, v,8), given the 
parameter f3\. But to solve these equations, it proves more convenient to regard (3\ as a 
variable, and 5 as a parameter; we adopt this point of view in what follows. 

There's another important structural aspect of equations (}28|) - (j3T| . namely, that (}30|) is 
distinguished in that it involves only two unknown quantities. It has the form u = f{u,5) 
and can be solved numerically for a given 5. When a solution exists, all other variables (ci, 
j3i, v) can be generated parametrically from the (u, 5) pair. Thus, the problem of solving 
equations (f2~% j) - (j3~Tj) reduces to a root-finding exercise in one dimension instead of four. 

Figure ITU1 plots the graph of the difference f(u, 5) — u for several values of 5. The zeros of 
this graph correspond to the solutions of (J30|) . and yield the desired (u(S), S) pairs. These 
are then substituted into the remaining equations to obtain ci(5), Pi(8), and v (5), from 
which various quantities of physical interest can be derived. 

A. Calculation of /drift 

For example, we can use the perturbative solution to find /drift, the fraction of the system 
that is drifting. It is most convenient to calculate this quantity first in terms of S, and then 
later re-express it in terms of the more natural control parameter 

To find the drifting oscillators, recall from Fig. that the cutoff between the locked 
portion and the drifting portion occurs at the crossover value x = x c where R(x c ) = |A|. 




(28) 
(29) 




(30) 
(31) 
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-0.05 



FIG. 10: Roots of Eq. (1M0D for various values of 5. Red indicates negative 5 and blue positive 5. 
For 5 < —0.028 there are no roots; for —0.028 < 5 < 0, two roots; for < 5 < g, one root; for 
h < 5 < 0.196, two roots; and for 5 > 0.196, no roots . 



Substituting (JT7j) for R 2 and equating this to A 2 , we obtain 



|c| + 2Re (ca*) cos x c + \a\ cos x c = A 



Plugging in the ansatz and keeping terms up to order e 2 , we find 



(32) 



1 + 2 Cl e + (c\ + 2c 2 )e 2 + 2ue 2 cos x c = 1 + 2A x e + (A 2 + 2A 2 )e 2 



(33) 



Finally, because of (|23jl . this simplifies to 



A, 



cos x r 



■ft 



c 2 _ 5 
tt 



(34) 



Since the spatial domain of the ring has length 2ir and 2x c is the length of the region occupied 
by drifting oscillators, the fraction of the chimera state corresponding to drifting oscillators 
is /^ft = x c /tt, and hence 



/drift - ~ 



TT 



COS 



-1 



s 



(35) 



FigurelTTlplots the numerically computed /(j r ift(^) against Pi(5). The curve has a turning 
point at (3\ ~ 0.22, when about 44% of the system is drifting. Presumably, this turning point 
stems from a saddle-node bifurcation in the underlying dynamics. In our simulations, we 
only see the upper branch of this curve, suggesting that this corresponds to the stable version 
of the chimera state. The reciprocal of the critical f3\ is about 4.5, which is the slope of the 
dashed line shown in Fig. |H1 in excellent agreement with the boundary of the wedge found 
numerically. 
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FIG. 11: Fraction of chimera state consisting of drifting oscillators as a function of (3\. Solid line 
indicates stable chimera, dotted line indicates unstable. The maximum (3\ determines the line 
bounding the wedge-shaped existence region in Fig. |HJ 




Homoclinic 



Locked Uniform Drift 

FIG. 12: Diagram of bifurcations giving rise to the chimera state in u — 5 plane. Insets show 
average frequency A versus x. Please see text for definitions of the perturbative variables 5 and u, 
and for a detailed explanation of this figure. 

B. Birth and Death of the Chimera State 

Although the parametric dependence of /^rift seems to be conveniently expressed with 
respect to /3\, that representation conceals a lot. Several dynamically distinct states of the 
system are invisible because they are all squeezed onto a single point (/3 1; /drift) = (0, 1), as 
we'll see below. It's much more revealing to use 5 instead of f3\. 

Therefore, we now examine the system in the set of coordinates shown in Fig. with 
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u plotted vertically and 5 horizontally. This picture is a compendium of all the stationary 
states of the system — the stable and unstable chimera, along with other states that we 
haven't mentioned yet, which we call uniform drift, modulated drift, and homoclinic locked 
states. The virtue of this representation is that it allows us to see each bifurcation that 
occurs as the chimera state comes into existence and later disappears. Beginning at the 
origin and moving counterclockwise around the kidney-bean shaped cycle, we have: 

1. Homoclinic locked state: u = 5 = 0. Here, all the oscillators are locked in phase, and 
hence frozen in the rotating frame. Accordingly, the average frequency A(x) of the 
oscillators vanishes for all x, as shown in the inset. But one can show that this state 
is not linearly stable. In fact, the exact, non-perturbative counterpart of this state is 
the in-phase locked state (|2~3|) at the critical parameter value (3 = 0, where this state 
undergoes a homoclinic saddle-node bifurcation. 

2. Spatially uniform drift: For u = and small 5 > 0, the system has a stationary state 
in which all the oscillators drift in a way that varies strongly in time but remains 
uniform in x. The order parameter R(x) is independent of x and close to 1, meaning 
that the oscillators are nearly in phase for nearly all of the time. An individual 
oscillator executes a jerky motion around its phase circle, lingering near 9 = and 
then whipping around the rest of the phase circle back to this point. The associated 
plot of A(x) is flat because of the uniformity in x. 

In terms of the perturbative variables used in Fig. this state appears on the line 
u = with 5 > 0. Then (|?TT|) shows that v = as well; hence a = 0, to O (e 2 ). So 
(j!7|) implies that R(x) must be real and constant and (fTSj) implies that Q(x) = 0. 
Equation (|2lJj) tells us that such a state is possible only if =0, which suggests that 
one can find an exact, non-perturbative version of the uniform drift state when (3 = 0. 
Indeed, seeking a solution of the self-consistency equation (fT^j) with j3 = 0, Q(x) = 0, 
and constant R(x), one finds R 2 = A — ^/ A 2 — R 2 since the kernel G is normalized. 
Hence, along this line (3 = 0, Q(x) = 0, and R(x) = R = y/2A — 1. 

3. Onset of spatial structure: At the lower right corner of the kidney bean, the function 
f(u, 5) — u (Fig. [TUj) becomes tangent at the origin, introducing a new branch of 
solutions with O = and (3 = but with the coherence R varying spatially. This is 
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the birth of spatial structure in the system. It happens for 8 = ~. The non-perturbative 
generalization of this result is A = 2/(2 + A). (See Appendix 1X1) 

4. Modulated drift: Along the first curved branch, all oscillators continue to drift, but 
now there is spatial structure in the R(x) curve, leading to a modulated pattern of 
average velocities (see inset). However, the average angle Q(x) is still identically 0. 

5. Chimera birth: At the point where u — 5, the first locked oscillators are born. For 
the first time, v and (3\ become nonzero (see Eq. (|2*9 ]) .([3T |) ). Until this point, all of 
the states have been confined to the vertical axis of Fig. |H1 now we finally we move 
off the wall. The curve of average velocities A(x) touches the x-axis at a single point. 
Meanwhile, the system develops spatial structure in its average phase: B(x) is no 
longer identically zero. 

This bifurcation can be shown to occur at 5 = 16/(97r 2 ), by evaluating the integral in 
(fHUj) with u = 5 ; also, see Appendix |B] for an exact calculation of the chimera state 
at birth. 

6. Stable chimera: Along the top of the kidney bean, the chimera state is dynamically 
stable. After its birth from the spatially modulated drift state, it gradually develops 
an increasing fraction of locked oscillators as we move in the counterclockwise direc- 
tion. Locked oscillators correspond to the zero part of the A(x) curve (they appear 
motionless because the reference frame was chosen to co-rotate with them). 

7. Saddle-node bifurcation: As we continue to move counterclockwise, the value of f3\ 
grows (along with the fraction of locked oscillators), and reaches its maximum at the 
point where the stable chimera ceases to exist. The disappearance is a result of a 
saddle-node bifurcation — a collision with an unstable chimera state — and occurs when 
about 44% of the system is drifting. 

8. Unstable chimera: Along the unstable dashed branch, the fraction of locked oscillators 
continues to grow. But the value of f3\ now begins to decrease, indicating a movement 
back towards the wall in Fig. |HJ The system returns to its original state when all 
oscillators become locked, with 5 = u = Pi = and R=l. 
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Taking a step back, we can see an interesting message of Fig. H21 The stable and unstable 
chimera states are continuously connected though the branches of drifting states, shown in 
solid red and dotted green lines. If we had used the /3 X representation instead (as in Fig.FTTj). 
both of these connecting branches would have shrivelled down to a point. The two kinds of 
chimera states would seem disconnected in a way that they really aren't. 

VII. DISCUSSION 

Our main result in this paper is an exact solution for the chimera state, for the special 
case of a cosine kernel. That solution also shed light on the bifurcations which create and 
destroy the chimera. 

In retrospect, it's not surprising that a cosine kernel would make the self-consistency 
equation (JTUJ) as tractable as possible, because the right hand side of (JTUJ) is a convolution 
integral, and trigonometric functions behave nicely under convolution. For this reason, it 
should be straightforward to extend the calculations to include more harmonics in G. Using 
the same argument as in Section IIVI one can see that the exact solution for the order 
parameter (JT1J) will have the same number of harmonics as G has. This approach would 
then give a systematic way to solve the self-consistency equation for any kernel representable 
as a finite Fourier cosine series. By taking more and more terms, this approach also gives a 
way to approximate results for any even kernel, as long as it is representable by a Fourier 
series. 

Unfortunately, the trick of choosing a special kernel may not work as well in two (or three) 
spatial dimensions. That could limit the applicability to two-dimensional chimeras, such as 
the novel spiral waves computed numerically in [10.]. Nevertheless, the idea of seeking a 
tractable kernel that can simplify the problem may itself be useful. 

Another caveat is that, despite its usefulness as a mathematical tool, the perturbative 
approach adopted here does not give a rigorous understanding of the bifurcations in the 
original problem. One would like to understand the bifurcation scenarios for all values of 
the coupling parameter A, which essentially measures how far the nonlocal coupling deviates 
from strictly global coupling. In Appendices El and El we show two results along these lines. 

One interesting aspect of the perturbative approach is that it draws our attention to the 
special parameter values A = 0, (3 = (or equivalently a = tt/2). Here the system has 
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global cosine coupling and is known to be completely integrable |23ll24j|. So in a sense, what 
we have done in this paper is perturb off this extremely degenerate system, which raises the 
question of whether other, unforeseen attractors might also lurk nearby, for different choices 
of initial conditions. 

The surprising nature of the chimera state makes us wonder if it could be created artifi- 
cially in a laboratory experiment, or possibly even occur naturally in some system. 

As a first attempt to judge whether this might be possible, we tried to integrate the 
phase equation with slightly nonuniform frequencies u>i, to mimic the inhomogeneities 
that would occur in any real system, and to test whether the chimera is an artifact of 
assuming identical oscillators. We added a uniform random variable r G [—B,B] to the 
native frequency u for each oscillator, and we found the chimera state persisted, as long as 
B was not too large — less than about 4% of A (the frequency difference between the locked 
oscillators' Q and the mean natural frequency u). This estimate should be conservative 
when compared with a presumably more realistic Gaussian random distribution of uj{. 

There are several possibilities for experimental systems where the distinctive effects of 
nonlocal coupling, including the chimera state, might be observed. Laser arrays seem to 
be goo d candidates. In some cases, such as semiconductor arrays with evanescent coupling 
iZQ] they are governed by equations similar to (JTJ), though these are usually approximated 
as nearest-neighbor. Likewise, phase equations of this form arise in the description of coupled 
electronic phase-locked loops, and superconducting arrays of Josephson junctions jsoL bd] . 
Finally, an idealized model of biochemical oscillators, coupled by a diffusible substance that 
they all produce, can give rise to an effectively nonlocal coupling and chimera states; indeed, 
this was the motivating example that led Kuramoto and his colleagues to their discovery. 

Whether or not the chimera state turns out to experimentally realizable, it is fascinating 
in its own right, as a strange new mechanism for pattern formation in spatially extended 
nonlinear systems. Its existence underscores how much still remains to be discovered, even in 
what would seem to be the simplest possible model of pattern formation: a one-dimensional 
collection of identical oscillators. 
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APPENDIX A: ONSET OF SPATIAL STRUCTURE 



We now show that the birth of spatial structure can be calculated non-perturbatively. We 
have already seen in Section IV 1 Bl that when (3 = 0, the system has an exact state of spatially 
uniform drift with constant coherence R(x) = V2A — 1 and average phase G(x) = 0. For 
this special state, the modulation amplitude a = and the mean level of the coherence 
c = R = \/2A — 1. At the bifurcation that creates spatial variation in the coherence, the 
real part of a becomes nonzero; at leading order in perturbation theory, this bifurcation 
takes place at 5 — |. Meanwhile, the imaginary part of a remains zero, which means that 
e(x) = still holds. 

To generalize this result to the non-perturbative case, we seek conditions for a second 
branch of solutions to bifurcate off the uniform drift state. Since A = (1 + c 2 )/2 for the drift 
state, we consider a slight perturbation 



A = (1 + c 2 )/2 + 7? 



(Al) 



where i] is a small deviation. Also, since a = for the uniform drift state, we may take a 
itself as a small deviation. Plugging all this into the self-consistency equation (}2*Tj) gives: 



'1 + c 2 )/2 + r] - [((1 + c 2 )/2 + T]) 2 - c 2 - 2accosa; - a 2 cos 2 a;] 



(A2) 



\ c + a cos x 

Now expand in a two- variable Taylor series for small i] and a, and integrate over x to obtain: 
2c 







-r} + 0[r} 



c(c 2 -3)(c 2 + 1) c(c 6 - 5c 4 + 19c 2 + 9) / 2 
H ~ —c V + 0[r] 



2(1 -c 2 ) 3 



(1 -c 2 ) 5 



Repeating the approach for the second self-consistency equation (|22|) gives: 
"2(c 2 - 1) + A(c 2 + 1) A(c 4 - 4c 2 - 1 







+ 



-r) + 0(r) 



+ 0(a 3 ) . 



(A3) 
(A4) 



2(c 2 - 1) (c 2 - l) 3 

To locate where another branch of solutions bifurcates off the uniform drift solution, we 
inspect the linearization of the algebraic system above, given by the Jacobian matrix 

2c 



(c 2 -l) 2 







2(c 2 -l)+A(c 2 +l) 
2(c 2 -l) 



(A5) 



If the determinant of the Jacobian is nonzero, the implicit function theorem tells us that 
no other solutions exist nearby. Hence, the existence of a continuously bifurcating branch 
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requires that the determinant vanish. Setting the determinant equal to zero yields c = 
or 2(c 2 — 1) + A(c 2 + 1) = 0. Plugging in the value of c about which we're linearizing, 
c = V2A — 1, and solving for A finally gives the bifurcation condition 

^~ 2 - (A6) 

To compare this with our earlier result from first-order perturbation theory, we express 
the perturbative variable 5 at this critical point by using its definition from ()27j) above and 
the property in (|23|) . Since A — c = 5e 2 = 5 A 2 (ignoring higher order terms), we have 
S = So 

Sr. = 




0(A A ) , (A7) 
which agrees with our perturbative prediction that S c = | in the limit that A — » 0. 
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APPENDIX B: BIRTH OF THE CHIMERA STATE 



When /3 = 0, it is possible to calculate the chimera state exactly, at the moment of 
its birth from a spatially modulated drift state. Recall that all states of pure drift satisfy 
Q(x) = 0, and equivalently, that a has zero imaginary part. Hence we can seek solutions of 
the algebraic self-consistency equations with real values of a and c, for (3 = 0. At the onset 
of the chimera, the first locked oscillators are born. As suggested by Fig. |3 this occurs when 
the graph of R(x) intersects the horizontal line R = A tangentially. 

Therefore the bifurcation condition is A = R max = c + a. Plugging this into (|T3jl and 
(jTfijl . and using Q(x) = and (3 = 0, we find that the self-consistency equations become 



' c + a — \/(c + a) 2 — (c + a cos x) 2 \ 

1 i L\ (Bl) 

c + a cos x / 



and 



/ c + a — J(c + a) 2 — (c + a cos x) 2 \ 

A{ ^ -cosx) . (B2) 

\ c + a cos x / 



Note that both of these expressions can be rewritten solely in terms of the ratio a/c, which 
suggests a neat way to solve them parametrically. 

Set s = a/c and substitute into (jBlj) above, which becomes 



1 + 8- J(l + S) 2 - (1 + SCOSX) 2 



\ 1 + s cos x I 

Ms) . (B3) 



So we can also write a = sc = sf\(s). 
Similarly, the a equation (|TB|) becomes: 

a = A(h cosx) 



s cosx) 2 



Jl + 8-y/(l + 8 y-{l 

= A { cos x 

\ 1 + s cos x I 

= Af 2 {s) . (B4) 

All other quantities of interest can also be expressed in terms of s. For instance, we can 
now substitute a = sc = sfi(s) into (jB4|) and solve for A(s) = a/fz(s) = sfi(s)/ f2(s). 
Likewise, A = c + a=(l + s)fi(s). In summary, the incipient chimera state can be written 
exactly in parametric form, as follows: 

c = fi(s) 
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a = sfi(s) 

f2[s) 

A = (1 + s)/i(s) . (B5) 

Since A is a control parameter of the original equations (the only free one after (3 has 
been chosen to be zero), it is desirable to reparametrize this solution in term of A. To do 
that, we invert A(s) in ()B5j) to obtain the following series expansion for s(A), 

16 a2 16 /3vr 2 -32\ 43 „ , , 4 

s ~ — -A — A 3 + O (A 4 

9tt 2 27 V 27vr 4 / V 

and use that to rewrite the newborn chimera in terms of A: 

16 , 8 /5tt 2 -32\ a2 . , 3 

c ~ 1 -A + - ) A + O (A 3 

3tt 2 9 1 vr 4 / V 



16 , 16 /3vr 2 - 16 



a ~ 



9tt 2 27 1 vr 4 / V 



16 . 8 /7tt 2 -32\ a2 _/ . 3 n 

Notice that this has exactly the form of the ansatz we postulated in ([24)1 . based on numerical 
experiments. As expected, it satisfies Ai = C\ as in (J25|) and gives A 2 — C2 = 16/(97r 2 ). 
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